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ABSTRACT 

A  survey  of  numerical  results  of  low-density  transitional  regime 
drag  coefficients  and  number  density  distributions  about  slender,  cold- 
wall,  conical  vehicles  at  hypersonic  speeds  is  presented.  The  transi¬ 
tional  regime  results  have  been  obtained  from  numerical  solutions  of 
Rosen's  restricted  variational  principle  representation  of  the  Boltzmann 
equation,  which  has  been  shown  to  be  equivalent  to  Galerkin's  method  of 
solution  of  operator  equations.  Transitional  regime  drag  coefficients 
are  compared  with  near  free  molecule  and  viscous  interaction  predictions 
of  drag  coefficients.  Number  density  profiles  for  the  transitional 
regime  are  also  given  and  compared  with  known  qualitative  features  of 
near  continuum  flow  fields  to  assess  the  validity  of  results  obtained  by 
the  variational  method. 

Results  indicate  that  the  drag  coefficient  varies  smoothly  from 
known  viscous  interaction  predictions  and  reaches  the  limiting  value  of 
2  when  the  freestream  Knudsen  number  is  on  the  order  of  unity  for  the 
class  of  slender  cones  considered.  A  decrease  of  cone  half-angle,  with 
Knudsen  number  fixed,  results  in  a  steeper  rise  of  the  drag  coefficient 
curve  toward  the  limiting  high  Knudsen  number  value  of  2.  For  cone  half- 
angles  on  the  order  of  5  degrees,  the  drag  coefficient  curve  exhibits  an 
overshoot  over  the  free  molecule  value. 
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I .  INTRODUCTION 

3 

In  Refs.  1  and  2,  Rosen's  restricted  variational  principle  method 

of  solving  the  Boltzmann  equation  was  used  to  determine  the  hypersonic 

low-density  transitional  regime  flow  over  conical  bodies.  Number  density 

profiles  as  well  as  drag  coefficients  were  calculated  for  Knudsen  numbers 

-2  1 

in  the  transitional  regime  between  10  and  about  10  for  the  case  of  an 

adiabatic  wall,  and  fhe  results  compared  with  those  yielded  by  inviscid 
4  \  5 

continuum  theory,  free  molecule  theory,  the  near-free  molecule  predic¬ 
tions  of  Laurmann,^  the  viscous  interaction  theory  of  Ellinwood  and 

7  8 

Mirels,  and  the  experimental  results  of  Hickman. 


The  present  note  presents  further  results  obtained  by  using  this 

variational  method  to  determine  number  density  profiles  and  dra  coeffi- 
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cients  in  the  Kn^  regime  between  10  and  10  ,  under  typical  reentry 
conditions  of  hypersonic  flow  over  cold-walled,  slender  cones.  Since 
Refs.  1  and  2  fully  present  the  theoretical  aspects  of  the  variational 
method,  only  numerical  results  are  given  in  this  note.  The  evolution  of 
number  density  profiles  with  decreasing  Knudsen  number  is  followed  down 

_3 

to  low  Knudsen  numbers  (Kn  =10  )  to  determine  whether  the  variational 

00 

method  is  capable  of  reproducing  certain  well-known  qualitative  details 
of  continuum  flow  fields.  Drag  coefficient  and  number  density  distribu¬ 
tion  results  are  given  for  5-  and  10-degree  cones  to  study  the  effects 
of  body  geometry  and  Knudsen  number  on  low-density  flow  fields,  and 
qualitative  trends  of  these  results  are  discussed. 


Abbreviated  Kn  and  defined  as  the  free  stream  mean  free  path  normalized 

00 

to  body  length. 
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II.  RESULTS 

Figure  1  represents  a  composite  of  predictions  for  the  theories 

given  in  Refs.  1-7  for  a  M  =  16  flow  over  a  cold  wall  10-degree  half- 
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angle  cone  in  the  Knudsen  number  range  between  10  and  10  (where  is 

the  free  stream  Mach  number) .  The  viscous  interaction  predictions  of 

Ellinwood  and  Mirels^  overlap  the  C  predictions  given  by  the  variational 

D  -2  -2 

method  for  Kn  in  the  range  between  10  and  3  x  10  .  Further,  the 

variational  method  results  qualitatively  agree  with  the  trend  of  Laurmann's 
near  free  molecule  C_.  variation  for  Kn  greater  than  1.  The  C_  variation 
predicted  by  the  variational  method  provides  a  smooth  transition  between 
the  limits  of  applicability  of  these  two  theories.  Figure  2,  which  gives 
a  summary  of  predictions  for  the  5-degree  half-angle  cone  case,  indi¬ 
cates  that  0^  increases  more  rapidly  for  the  5-degree  cone  for  Kn^  greater 
than  about  10  ^  than  for  the  10-degree  cone.  Further,  a  small  over¬ 
shoot  over  the  free  molecular  value  is  evident  both  from  the  near  free 
molecule  and  variational  theories,  although  the  overshoot  is  less  than 
might  be  expected  by  the  trend  of  the  near  free  molecular  theory.  Toward 
the  lower  Knudsen  number  end  of  the  curve,  the  region  of  overlap  of 

variational  and  viscous  interaction  theories  appears  to  begin  for  Kn 
_2 

below  about  10  .  For  both  the  5-  and  10-degree  cones,  is  close  to 

the  free  molecular  limiting  value  for  Kn^  of  about  1. 

Figures  3,  4,  and  5  represent  the  evolution  of  number  density 

profiles  with  increasing  Knudsen  number  for  the  10-degree  cone  and  give 

-2  -1 

number  density  profiles  for  Kn^  =10  ,10  ,  and  1,  respectively.  As 

Kn^  decreases  the  disturbance  region  appears  to  narrow  (compare  Fig.  3 

to  Fig.  5)  in  the  vertical  direction.  For  higher  Kn^,  for  instance 

Kn  greater  than  10  \  the  sensible  disturbance  region  appears  confined 

* 

to  distances  less  than  a  free  stream  mean  free  path;  only  for  Kn^  below 
10  ^  do  the  disturbance  regions  appear  to  exceed  one  free  stream  mean 


* 

The  mean  free  path  referred  to  here  is  based  on  free-stream  conditions; 
within  the  disturbed  flow  region  the  mean  free  path  is  less. 
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MOLECULE  LIMIT 


Figure  1.  Drag  Coefficient  as  a  Function  of  Kn^,  =  10  degrees 
M  =16,  cold  wall  conditions 
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Figure  3.  Number  Density  Profiles  =  16,  =  10 
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Figure  4.  Number  Density  Profiles 
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Figure  5.  Number  Density  Profiles  M  =  16,  iji,  =  10°,  Kn  =1.0 
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free  path  in  vertical  extent.  The  appearance  of  a  number  density  maximum 

in  the  flow  adjacent  to  the  cone  is  apparent.  For  Kn^  equal  to  10  1  and 

10  2  the  maxima  are  at  roughly  the  same  height  above  the  cone  surface, 

_2 

although  for  Kn^  equal  to  10  the  increasing  importance  of  gas  particle 
collisions  appears  to  cause  a  decrease  in  the  particle  population. 
Although  the  vertical  height  of  the  number  density  maxima  appears  to  be 
the  same,  the  number  of  mean  free  paths  away  from  the  cone  surface  of 
these  maxima  continually  increases  as  Kn^  decreases.  In  all  cases,  the 
appearance  of  number  density  maximums  occurs  towards  the  rear  of  the  cone. 

_3 

Figures  6,  7,  8,  and  9  are  number  density  profiles  for  Kn  =  10  , 

-2  2  “ 

10  ,  10  and  10  ,  respectively,  and  represent  the  evolution  of  number 

density  profiles  for  a  5-degree  half-angle  cone.  For  Fig.  6,  the  Knudsen 

_3 

number  is  10  ,  which  is  close  enough  to  continuum  flow  conditions 

(Kn^  ->•  0)  that  qualitative  features  of  the  continuum  shock  layer  should 
be  apparent  if  the  variational  theory  is  valid  to  such  low  Kn  values. 

OO 
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Figure  7.  Number  Density  Profiles  ■  16,  <|j  *=  56,  Kn^  «=  0.01 
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The  results  show  that  the  diffusive  number  density  profiles  characteristic 

of  higher  Kn^  values  (Fig.  7  for  example)  develop  an  increasingly  sharp 

gradient  with  decreasing  Kn  (Fig.  6),  typical  of  a  fully  developed  shock 

-3 

layer.  The  amplitude  of  the  number  density  maximum  in  the  Kn^  =  10 

case  (Fig.  6)  is  ri/ri  =  2.12;  ThOs  compares  favorably  to  the  post-shock 

number  density  predicted  from  classical  inviscid  continuum  theory  of 

p/n  =  2.6.  Further,  a  line  drawn  through  the  number  density  maxima 

compares  favorably  to  the  inviscid  continuum  shock  position  (shown  in 

Fig.  6).  For  both  the  5-  and  10-degree  cone  cases,  a  definite  evolution 

occurs  from  the  diffusive  small  gradient  profiles  characteristic  of  high 

Kn^.  Comparison  of  Fig.  7  to  Fig.  3  indicates  that  the  higher  cone  angle 

yields  higher  maximum  number  density  values  in  a  narrower  disturbance 

region  than  for  the  5-degree  cone.  The  latter  case  is  characterized  by 

more  diffusive-type  profiles  than  the  10-degree  cone  but  still  exhibits 

the  qualitative  feature  of  having  number  density  maxima  away  from  the 

wall.  Towards  the  tip  region,  density  gradients  are  small;  larger 

gradients  are  only  developed  towards  the  rear  portions  of  the  cone  for 
_2 

Kn^  below  about  10  for  the  5-degree  cone  and  below  about  2  for  the 
larger  10-degree  cone.  In  all  the  cases  calculated,  wall  number  density 
exceeds  the  free  stream  value.  This  is  due  both  to  the  fluid  compression 
effect  which  leads  to  an  increase  in  wall  density  with  increasing  cone 
angle  (compare  Fig.  3  to  Fig.  7  for  instance),  and  to  the  effect  of  the 
cold  wall,  which  becomes  more  efficient  in  raising  the  wall  number 
density  level  as  Kn^  decreases  (as  wall  velocity  slip  decreases  as  Kn^ 
decreases;  compare  Fig.  9  to  Fig.  6,  for  example). 

Comparison  of  disturbance  regions  associated  with  the  high  Kn^  flow 
over  slender  cones  (compare  Figs.  6  to  9  for  the  5-degree  cone  to 
Figs.  3  to  5  for  the  10-degree  cone)  shows  that  larger  values  of  number 
density  exist  further  from  the  cone  surface  for  the  more  slender  cones. 
This  effectively  results  in  an  increased  probability  of  collisions  of 
free  stream  particles  with  the  larger  cloud  of  particles  surrounding  the 
more  slender  cone.  The  occurrence  of  large  number  density  regions  away 
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from  a  body  surface  and  the  consequent  increase  of  scattered  particles 

9 

to  the  body  surface  for  high  Kn^  flows  have  been  noted  by  Bird  to  result 
in  a  near-free  molecule  CD  overshoot  above  the  free  molecule  limiting 
value.  The  presence  of  the  overshoot  for  the  5-degree  cone  is  consis¬ 
tent  with  the  larger  disturbance  region  for  this  geometry  compared  to 
the  10-degree  cone  case,  for  which  no  overshoot  exists. 

A  further  comparison  of  the  free-stream  Mach  number  effect  on  the 
extent  of  the  number-density  disturbance  region  is  obtained  by  comparison 
of  Fig.  10  to  Figs.  4  and  5.  In  Fig.  10  the  Mach  number  is  5.4;  in 
Figs.  4  and  5  it  is  16.  For  the  lower  Mach  numbers,  the  vertical  extent 
of  the  disturbance  region  is  on  the  order  of  10  to  15  times  that  for  the 
higher  Mach  number.  In  addition,  the  profiles  of  Fig.  10  exhibit  a  more 
diffuse  character  due  to  the  greater  importance  of  thermal  motions  of 
gas  particles  at  lower  Mach  numbers. 

The  survey  of  number  density  distributions  and  drag  coefficients 
given  in  this  report  is  meant  to  provide  details  of  transitional  regime 
flow  over  conical  vehicles  at  typical  reentry  conditions.  Little  theo¬ 
retical  or  experimental  work  has  yet  been  done  on  transitional  regime 
flow  fields  over  cones,  although  viscous  interaction  and  free  molecule 
theories,  which  give  asymptotic  limits  that  transitional  regime  results 
must  agree  with,  are  available.  In  comparison  with  these  theories,  the 
variational  method  yields  smooth  transitions  of  drag  coefficient  vari¬ 
ation  between  known  correct  limiting  cases;  additionally,  details  of 
continuum  and  free 
variational  method 


molecule  flow  fields  are  likewise  obtained  from  the 

-3  2 

for  cases  of  Kn  =10  and  Kn  =  10  ,  respectively. 
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